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Abstract 



The exact solution of a general Z(4) gauge Potts model with a single and 
double plaquette representation of the action is found on a subspace of gauge- 
coupling parameters on the square and triangular lattices. The two Ising-type 
critical lines of a second-order phase transition for the model on a square lat- 
tice are found. For the model on a triangular lattice the two critical surfaces 
of an Ising-type and two nontrivial lines of a second-order phase transition 
with different critical behavior than on the critical surfaces are found. It is 
shown that a two dimensional (2D) general Z(4) gauge Potts model with sin- 
gle and double plaquette representation of the action and a 2D spin-| Ising 
model belong to the same universality class. 
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I. INTRODUCTION 



Confinement-deconfinement phase transition in gauge theory is one of the challenging 
problems in modern physics. Lattice-gauge theory is a gauge- invariant nonperturbative tool 
of regularization of gauge action to avoid the ultraviolet divergence in the theory. Although 
this regularization procedure is not unique, different ways of defining gauge theories on the 
lattice should lead to the same physics in each case when the continuum limit of vanishing 
lattice spacing is taken. From a theoretical point of view, investigations with different lattice 
actions will enable a deeper understanding of the physics of confinement and other related 
problems in QCD. The first lattice-gauge model with a single plaquette representation of 
the action has been introduced by Wilson, 



s = -pJ2 ReU i 



p ■ 



where U p denotes the usual plaquette variable, the product of link gauge fields around a 
plaquette It was expected that non-Abelian gauge theories, in general, do not have any 
phase transitions separating strong- and weak-coupling regimes. Therefore, confinement, 
explicitly shown on the lattice in the strong-coupling region, should persist also in the con- 
tinuum limit. Later, Bhanot and Creutz [[| extended the form of the Wilson action by 
adding an adjoint coupling term. Using Monte Carlo simulations it was shown that confine- 
ment could survive even through the phase diagram of the mixed action and the so-called 
bulk (volume) phase transitions separating strong- and weak-coupling regions exist ||@]. 

The choice of action is still far from unique. Recently, several improved actions have been 
proposed as a way of reducing scaling violation in the approach to the continuum limit from 
a lattice action. Among them the Symanzik-Wiesz action constructed from a combination 
of (1 x 1) and (1 x 2) Wilson loops, the Bhanot-Creutz action 0j, the tadpole-improved 
actions and the g-state gauge Potts model with a single and double plaquette form of 
action 0. 

The Monte Carlo analysis of the g-state Potts model with a single and double plaquette 
form of action |J showed that for d = 3, q — 2, first- and second-order transition lines; and 
for d = 2, q = 3, 4, the second-order; and for q = 5, first-order transitions are existed, which 
is in good agreement with the analytical results for d = 2. 

The lattice-gauge models with double plaquette interaction terms in the action were 
proposed and studied in three dimension and four dimension by Edgar and Bhanot et 
al. ||. Turban investigated the two dimensional 2D gauge model with the global Z{2) 
symmetry on a rectangular lattice ||. He reduced it to the usual spin-| Ising model on 
a square lattice and obtained a point of a second-order phase transition. The Z(3) gauge 
model on the flat triangular and square lattices with double plaquette representation of the 
action was investigated by Ananikian and Shcherbakov [|10j. It was reduced to the spin-1 



Blume-Emery-Griffiths (BEG) model [11] and an Ising-type critical line of a second-order 
phase transition was found on a subspace of the interaction constants [ |T2| . 

The fact that lattice gauge theories could be mapped to the classical spin systems is well 
known. For example, Wilczek and Rajagopal |13[ showed that in QCD with two flavors of 



massless quarks, the chiral phase transition is in the same universality class as the classical 
0(4) Heisenberg antiferromagnet and they also established a dictionary between QCD and 
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the magnetic system (see also Ref. |14]]). Okawa JT5], using Monte Carlo renormalizat ion- 
group methods, showed that a (3 + l)-dimensional SU(2) lattice-gauge theory and a three- 
dimensional Ising model belong to the same universality class. 

In this paper we consider the generalized Z(4) gauge Potts model with a single and double 
plaquette representation of the action on the square and triangular lattices. We found an 
exact analytical solution of this model on a subspace of gauge-coupling parameters. Using 
duality transformation |Ll| and exact results for the spin-| Ising model on the square and 



honeycomb lattices [17].[T8|], we investigated the critical properties of the gauge theory. We 



showed that a two-dimensional generalized Z(4) gauge Potts model with single and double 
plaquette representation of the action and a two-dimensional spin-| Ising model belong to 
the same universality class. 

The paper is organized as follows. In Sec. II we define the model under consideration 
and present obtained results. Section III contains some concluding remarks. 



II. THE Z(4) GAUGE POTTS MODEL 

The most general form of the Z(4) gauge Potts model with a single and double plaquette 
representation of the action is defined through the action 

S Gauge = ~ ^ X) A/* 8u Pi >V 8u Pj ^"EE^^P,^ ' ( 1 ) 

(p iPj ) V,z Pi z 

where the outer summation in the first term runs over all nearest-neighbor plaquettes and 
in the second one is over all plaquettes of the lattice. The indexes y and z of the inner 
summations run over the group Z(4). The U p = TibedpUb denotes the ordered product of 
link gauge fields U^s around an elementary plaquette. Each link variable Ub takes the value 
exp(ikir/2) G Z(4), k = 0,1,2,3. 5 is the standard Kronecker symbol and f3 yz ,f3 z , y,z G 
Z(A) are coupling parameters. 

From the obvious identity for the Kronecker symbols 

Sv Pi ,i + du Pi , Zl + 8u Pi , Z2 + fiu Pi ,z3 = 1 > 

where G Z(A), k = 1, 2, 3, and an assumption that the coupling parameters are symmetric 
under the transposition of the indexes y and z, we can reduce the number of independent 
coupling parameters and rewrite gauge action (P in the following form 

Scauge = ~ X [^11^,1^,1 + P22^U Pi ,z^U Pj ,Zi + Pm&U Pi ,z 2 ^U Pj ,z 2 + 
(PiPj) 

P12 (S UptA 5 Up . )Z1 + S Up ^ Zl 5 UpjA ) + f3 13 (d Upiil 5 Up . tZ2 + 5 UpijZ2 5u Pj ,i) + (2) 

U Pi ,zi&Vp.,Z2 + fiu Pi ,Z2$U Pj ,zi)\ + X iP^u Pi ,i + {h8u Pi ,zi + Pz&u Pi ,z^) ■ 

Pi 

The partition function for this gauge model is defined as a sum of Boltzmann weights 
exp(S Gauge) over all configurations of the gauge variables {U}, 
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Z Gauge — E eX P [ _ $ Gauge] ■ 
{U} 



(3) 



To establish the connection between this gauge model and the spin-| Ising model, we intro- 
duce spin variables Si in the sites of the dual lattice such that 

Si = 2(^,20 - &u PvZl ) + l(8u Pi ,z 2 - $u Pi ,z 3 ) , 

Sf = i(S Upi>Z0 + 8 Upi>zl ) + l(6 Upitta + S UpvZ3 ) , (4) 
Using this substitution, we can rewrite gauge action (fj) in terms of new spin variables 



Ssni 



E 

(ij) 



J SiSj + K S*S] + L SfS* + y (SiS] + SfSj)+ 



J2(hSi-ASf + h 3 Sl 



(5) 



where constants J, K, L, M, Mi, M2, h, A and /13 are linear combinations of the original 
gauge-coupling parameters: 



J 
K 
L 

M 
Mi 
M 2 
h 
A 
h 3 



3 81 3 

32^ 13 + 64^ 22 + 32 ^ 23 



576 







o /?13 + 777/322 

8 lb 



7P23 + 77/^33, 
8 lo 



1 1 

— 011 + -012 - 

1011 — 7,011 ~ 7:013 + ~7022 + 7:023 + 777^33; 
4 2 O 4 3o 

9 ~ 9 a ^ o o 1 

— gPll + 7^12 + 7^13 — gP22 — J^023 ~ 77^P 33 , 

9 7 913 1 

"IT 11 + i/ 13 + W 22 " 24^ 23 " IT 33 ' 

7Pll — 7Tpl3 — -7022 + -023 + 77^33; 
4 3 4 Iz 

81 o 3 ^ 81 o 3 a 1 a 
011 013 022 H #23 H #33 

128 32 M 3 128 64 M 3 384 P33 



33 • 



+ 8 A " 8 A " 24*' 



^011 + 7^012 ~ 7^013 + 1^022 ~ 7^023 + 7^033 



+ \01 + \02 ~ \fc 



9 19 11 

7 ( ~32 Pn + 8 As + 32^ + l/ 23 ~ 96^ 33 



\01 + \02 + \03. 



where 7 is the coordination number of the dual lattice. Action (^) coincides with the 
Hamiltonian multiplied by 1/ksT of the spin-| Ising model [|T7 JT8 ] . Thus, the partition 
function of generalized gauge Potts model (§) is equal up to the factor to the partition 
function of the spin-| Ising model defined on the dual lattice 



iy / \a{"f)N yDual 

A Gauge ~ ^ ^ Spin > 



(6) 
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where 



yDual 
^Spin 



2^exp 

{s} 



■S. 



Spin] 



and 0(7) is a constant that depends on the coordination number of the lattice, a = 1 for the 
square lattice and a = | for the honeycomb one respectively. 

A factor A a ^ N has been included in Eq. (^) to take into account the difference between 
the number of gauge {U} and spin {S} configurations. To obtain the phase structure of 
this gauge model we will restrict ourselves to the spin-| Ising model on the square and 
honeycomb lattices. For coincidence with the spin-| Ising model the coefficients h, h 3 , Mi 
and M 2 will be set to zero. Thus, we can express the parameters of the spin-| model through 
the rest of the gauge couplings as follows: 



J 



K 



M 



A 



337 o 81 a 

P11 0i2 H 

288 32 H 288 



1 49 

KL3 + 7^7 P22, 

6b 



1 1 

o/3ll + 0A2 



'-A 



13, 



1114 
25 9 1 14 

— T^PH + -P12 — — P13 — P22, 

do 4 do 9 
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9 9 
32 Al + 32 A2 



32 



13 



The general spin-| Ising model on a honeycomb lattice was investigated by Izmailian 



and Ananikian [18]. The model is described by the Hamiltonian 



t/ SiSj -\- k s\ s*j -|- sf Sj -\- (SiSj + sfs. 
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where (3 = l/k B T, Si = ±|, ±| is the spin variable at site i, and (ij) indicates the 
summation over all nearest-neighbor pairs of sites. Under the conditions 



tanh 2 (Ji) 
exp (-AK) 



tanh(J 2 ) tanh(Jo), 
cosh ( J 2 - J ) , 



(7) 



that in terms of gauge-coupling parameters are 



Jo — i/3n — 5A2 + 4/^22, 

Jl = 4 All — 4/^22, 

^2 = — 2 All + 9A13 + l/^22 



(8) 



the model transforms to the spin-| Ising model on the same lattice. The free energy and 
critical point for the spin-| Ising model on the honeycomb lattice in the limit of an infinite 

lattice are well known |l9| . Thus, using this result one can obtain the important thermody- 
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namic properties of the spin-| Ising model with Z(2) symmetry ||17fl . After substitution of 
J, K and L from Eq, 
spin- § 



into Eq. (ffl) we obtain the subspace in which the corresponding 



3 Ising model can be solved exactly, 
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where 



exp(-4i?) 




cosh4 - (/3i2 + /?i3 - /3u) , 



in the space spanned by /3i 2 , /?i3, and (3\. It is easy to see that the A surface of the critical 
points in Eq. (|K]) is defined only in the two regions of the (/3i 2 , ,$13) plane, 

(i) < tanh 1(022 - 2^2 + )9u) < 1/V3 and 
1A/3 < tanh | (/3 22 + 2/3 13 - /3 n ) < 1, 

(ii) < tanh |(/3 22 + 2/3 13 - /?u) < 1/V3 and 
l/v/3 < tanh | (/3 22 - 2/3 12 + /3 n ) < 1. 

For each set of f3 12 and /?i 3 , Eq. (p^0|) determines the unique value of A, except for the 
intersecting point of the two regions (i) and (ii) for which (3\ is an arbitrary. Thus, the A 
surface in Eq. (|ToD contains two nontrivial A-lines of critical points given by 

(a) fti = /3i3 = 2 ln(2 + V3), (3 12 = P22 = 0, ^-arbitrary, 

(b) Ai = P12 = As = 0, (3 22 = 2 ln(2 + v^), ^-arbitrary 

As shown in |18|] on the A lines the model exhibits a critical behavior different from the 
critical behavior elsewhere on the A surface. The phase transition is not marked with the 
logarithmic divergence of the derivative of the order parameter P, where P is the quadrupolar 
moment defined as 



This phase transition is associated with the logarithmic divergence in the specific heat. 

The area in the plane of coupling parameters exp(/5i 2 ) and exp(/3i 3 ), where the A surface 
exists is shown in Fig. 1. In the Appendix it is proved that this area is connected and there 
is no phase transition for T —>■ 00 and T —>■ 0. Hence, for all possible values of the coupling 
parameters there is only one finite critical value of the external field for which the phase 
transition is of the second-order except points (a) and (b) in the above equation. 
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The spin- 1 Ising model with Z(2) symmetry was investigated on a square lattice by 
Izmailian [|l7j. It was shown that this model is reducible to an eight- vertex model on a 
surface in the parameter space spanned by the coupling constants J, K, L, and M. It was 
also shown that this model is equivalent to an exactly solvable free fermion model along two 
lines in the parameter space. The two A lines of a second-order phase transition was found 
exactly in this model. 

In terms of our gauge theory, these A lines have the following form: 

(c) fri = Pu = Pn = 0, p 22 = 2 ln(l + y/2), ^-arbitrary, 

(d) p u = P l3 = 2 ln(l + V2), p 12 = p 22 = 0, ^-arbitrary. 

On these A lines our gauge theory exhibits an Ising-type second-order phase transition 
(logarithmic specific-heat singularity). 

Thus, we showed that there exists an area (lines for square lattice) where the 2D gener- 
alized Z(4) Potts gauge model mapped to the corresponding 2D spin-| Ising model. Hence, 
because of universality of critical indexes it follows that these two models have the same 
critical indexes and belong to the same universality class. 

III. CONCLUDING REMARKS 

In summary, we have found an exact analytical solution of the Z(4) gauge-lattice model 
with a single and double plaquette representation of the action by mapping it to the dual 
spin- 1 Ising model with Z{2) symmetry. For the model on the square lattice we found the 
A lines of the second-order phase transition with logarithmic specific-heat singularity. For 
the model on the triangular lattice we derived the two A surfaces of a second-order phase 
transition with a usual Ising-type singularity of the order parameter and two nontrivial 
A lines of critical points on which our model exhibits the critical behavior unlike critical 
behavior elsewhere on the A surfaces. We demonstrated that the 2D general Potts gauge 
model belongs to the same universality class as the 2D spin-| Ising model. 
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APPENDIX 

Here we present the analytic investigation of the A surface in the plane of the coupling 
parameters exp(/3i 2 ) and exp(/? 13 ) for triangular lattice. Let us make the following denota- 
tions x = exp(/3i 2 ), y = exp(/3i 3 ) and z = exp(p 22 ). After elimination of Pn from Eq. @ 
one obtains the first-order polynomial for z, 

z(l-2y + xy) + (2-x)(y-x) = 0. (Al) 

In terms of the variables x, y and z the conditions (i) and (ii) take the following form: 
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A(i) 1 < zy{2 - x)/x 2 < c 2 and zy/(2 - x) > c 2 , 
A(ii) 1 < zy/(2 - x) < c 2 and zy{2 - x)/x 2 > c 2 , 

where c = 2 + a/3. From A(i) and A(ii) it is easy to show that x < 2 for any values of z 
and y. Taking into account that for T — > x,y,z can take only values {0, l,oo} and for 
T^oo £ = y = z — 1, one can show that neither A(i) nor A(ii) are satisfied, hence there 
is no phase transition for T — > and T — > oo. Using this fact one can construct the area 
in the plane of coupling parameters exp(/3i2) and exp(/?i3) in which Eq. (Al) and one of the 
conditions (Ai) or (Aii) are satisfied. This area is shown in Fig. 1. 
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FIGURE CAPTION 



Fig 1. The area in the plane of the coupling parameters exp /? 12 and exp /? 13 where the critical 
A surface exists, i.e., conditions (H) and (i) or (ii) are satisfied. Points A and B correspond 
to projections of two nontrivial A lines of (a) critical points and (b) for which (3\ is arbitrary. 
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